Formation of atomic tritium clusters and condensates 
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We present an extensive study of the static and dynamic 
properties of systems of spin-polarized tritium atoms. In par- 
ticular, we calculate the two-body \F, tuf) = |0,0) s-wave 
scattering length and show that it can be manipulated via a 
Feshbach resonance at a field strength of about 870 G. Such 
a resonance might be exploited to make and control a Bose- 
Einstein condensate of tritium in the |0, 0} state. It is further 
shown that the quartet tritium trimer is the only bound hy- 
drogen isotope and that its single vibrational bound state is a 
Borromean state. The ground state properties of larger spin- 
polarized tritium clusters are also presented and compared 
with those of helium clusters. 
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In 1976 0, Stwalley and Nosanow suggested, based on 
statistical arguments, that the system of spin-polarized 
bosonic tritium atoms behaves "very much like 4 He" . To 
the best of our knowledge, their arguments have not yet 
been tested by a microscopic quantum mechanical treat- 
ment. A detailed theoretical study of spin-polarized tri- 
tium systems, namely spin-polarized atomic tritium clus- 
ters and optically-pumped tritium condensates, is the ob- 
jective of the present work. In this Letter, we present 
results for both structural and scattering properties of 
tritium dimers, trimers, and clusters. 

Pioneering experimental studies of the lowest quartet 
state of spin-polarized atomic trimers have been pursued 
recently for sodium and potassium Study of these 
trimers, which were prepared on the surface of large 4 He 
clusters, revealed that three-body effects are surprisingly 
important 0]. We are not aware, though, of any ex- 
perimental or theoretical studies of larger spin-polarized 
atomic cluster systems. Bosonic helium systems — i.e. 
liquid bulk 4 He,two-dimensional 4 He films, and finite size 
4 Hejv clusters j| — have, of course, been studied exten- 
sively. This Letter thus presents the first predictions for 
spin-polarized atomic clusters. In particular, we charac- 
terize spin-polarized tritium clusters [in the following de- 
noted by (Tf)jv] with up to N = 40 tritium atoms, and 
compare their energetic and structural properties with 
those of bosonic 4 He at clusters. We hope that this study 
will stimulate further experimental work. Of particular 
interest is the lowest quartet state of the tritium trimer, 
which we predict to be a Borromean or halo state, and 
tritium cluster formation in the presence of an external 
magnetic field. 

We also point out the possibility for creating an 



optically-pumped gaseous tritium condensate. Bose- 
Einstein condensates (BECs) are, to first order, well 
characterized by the two-body s-wave scattering length 
between two atoms. It will be shown below that the 
triplet two-body s-wave scattering length at of two tri- 
tium atoms is large and negative, implying an unstable 
condensate of spin-polarized tritium atoms. We find, 
though, that there is an unusually broad Feshbach reso- 
nance H for two high-field-seeking, \F, rap) — 1 0, 0) tri- 
tium atoms (F denotes the total angular momentum, and 
to p the magnetic quantum number of this state). 

While condensation of atomic hydrogen was realized 
experimentally in 1998 |6), it was a difficult experiment 
- at least in part because of hydrogen's small triplet 
scattering length that limits the utility of evaporative 
cooling. Nevertheless, owing to hydrogen's simplicity, 
it remains an important species to study. For instance, 
properties such as the interatomic potential and spin re- 
laxation rates can be obtained theoretically from first 
principles. Unfortunately, we find no Feshbach resonance 
at reasonable field strengths for hydrogen. Thus, the res- 
onance for tritium may permit faster condensation of a 
hydrogen-like atom, and allow for the formation of a sta- 
ble BEC of tritium atoms with controllable properties in 
an optical dipole trap Q. Formation of such a tritium 
condensate should enhance the lively interplay between 
theory and experiment. 

The behaviour of atomic tritium clusters and conden- 
sates is primarily determined by the two-body interaction 
potential for two tritium atoms, which is identical to that 
for two H or D atoms, except for the isotope-dependent 
adiabatic correction. Since there are only two electrons, 
these dimers are among the few for which highly accu- 
rate ab initio potentials are available. In the following, we 
concentrate on the singlet ground state (S = 0, X 4 E+), 
and on the triplet ground state (S — 1, fo 3 £+) of the 
tritium dimer. 

To construct the two-body S = and S = 1 Born- 
Oppenheimer interaction potentials for hydrogen and its 
heavier isotopes, highly accurate ab initio data for the 
short range part (|] that incorporate the mass-dependent 
adiabatic correction are connected smoothly with an 
analytical expression describing the long-range behav- 
ior [p|~pd| . This procedure results in six potential curves 
describing H2, D2, and T2 in their S — and S = 1 
states, respectively, which are then used in the radial 
Schrodinger equation describing the relative motion of a 
particle with reduced mass to/2. The mass dependence 



1 



of the singlet and triplet scattering lengths has been dis- 

13]. Here, we esti- 



cussed extensively in the literature 1 12 
mate the uncertainty of our two-body scattering observ- 
ables by solving the radial Schrodinger equation using 
both the reduced atomic mass and the reduced nuclear 
mass for each two-body potential described above. 
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FIG. 1. Tritium dimer triplet 6 3 Ej potential (solid line) 
together with the He dimer potential (dotted line) as a func- 
tion of the interparticle distance r. The inset compares the 
tritium dimer triplet potential (S — 1, solid line) with the 
tritium dimer singlet ground state potential (S = 0, dashed 
line) . Note the different vertical scales of the main figure and 
the inset. 

Figure [j] compares the (Tf)2 potential (solid line) with 
the similarly shallow He2 ground state potential, LM2M2 
from Aziz and Slaman jL4j (dotted line). The Tf dimer 
and the He dimer potentials have a well depth of D e = 
— 4.6cm -1 and — 7.6cm -1 , respectively. Note that on the 
scale shown in Fig. [j} the Hf and Df dimer potentials 
would be indistinguishable from the T| dimer potential. 
The minimum of the (T| )2 potential lies at a significantly 
larger interparticle distance (r e = 7.8a.u.) than for the 
He dimer (r e = 5.6a.u.). Given the lighter mass of the 
tritium atom (to = 5496. 9a. u.) compared to that of the 
4 He atom (to = 7296. 3a. u.), it is not surprising that the 
Tf dimer is not bound, even though the tritium van der 
Waals coefficient Cg = 6. 499a. u. is larger than that for 
He, Cq — 1.367a.u. (recall that the 4 He dimer binding 
energy is only —9.1 x 10 -4 cm -1 For comparison, 

the inset of Fig. |l| shows the tritium triplet potential 
(solid line) together with the tritium singlet potential 
(dashed line). The singlet curve is almost four orders of 
magnitude deeper than the triplet curve and supports 27 
vibrational s-wave bound states. 

In agreement with values tabulated in the litera- 
ture |n|, we calculate the two-body s-wave triplet scat- 
tering lengths a t for (H|)2 to be a t — 1.33a.u. [1.33a.u.], 
and for (Df) 2 to be — 6.89a. u. [— 6.88a. u.], using the re- 
duced atomic [reduced nuclear] mass. From symmetry 
considerations, the s-wave scattering length for (D|)2 is 
not an observable, and is given here for diagnostic pur- 
poses only. For tritium, we predict a positive singlet scat- 
tering length, a s = 34.6a.u. [35.8a.u.], and a large nega- 



tive s-wave scattering length, a t = — 82.1a. u. [— 81.9a. u.]. 
The Tf dimer does not possess a bound state, but this 
large negative at indicates that it is only "slightly short 
of binding" . 

Despite the fact that the tritium triplet scattering 
length is negative — implying an unstable condensate 
— it may be possible to form a stable tritium condensate 
utilizing a Feshbach resonance. Coupled-channel scatter- 
ing calculations that couple the singlet and triplet sub- 
spaces reveal such a Feshbach resonance, i.e. a diverging 
scattering length for two atoms characterized by quan- 
tum numbers F and tuf, as a function of the magnetic 
field strength. The coupling arises through the atomic 
hyperfine interaction, and has to be accounted for by 
an effective two-atom Hamiltonian Jn|. Feshbach reso- 
nances have been observed experimentally for 23 Na and 
85 Rb H among others. The latter paved the way for a 
series of important BEC experiments, some of which en- 
ter the large interaction regime fl7j] while others probe 
the collapse regime O]. 
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FIG. 2. a(00 + 00) scattering length (diamonds; using the 
reduced atomic mass in the coupled-channel calculation, see 
text) as a function of the magnetic field strength B (dotted 
lines are shown to guide the eye). The solid line describes the 
behavior for the range B £ [400, 1300G] well using the fol- 
lowing parameterization, a(00, 00) = o.bg(1 — &/(B — Br)) 
with cleg = -37.9a.u., A = -1238G, and Br = 870.8G; 
however, the fit is inaccurate at low fields. Inset: Threshold 
energies in GHz as a function of magnetic field B in Gauss. 
The assignment of quantum numbers is approximate, except 
for 1 1 , 1 ) + 1 1 , — 1 ) , which is an exact eigenstate without ad- 
mixtures. 

We find that collisions between two tritium atoms, each 
in their \F, rap) — 1 0, 0) state, results in a scattering 
lenth of a(00,00) = 57a. u. [64a. u.] for zero magnetic 
field strength B, using the atomic [nuclear] mass. As 
B increases, the scattering length a(00 + 00) rises and 
eventually goes through infinity across a broad Feshbach 
resonance centered at B = 870G [810G] (see Fig. |). We 
also looked for, but could not find, a similar resonance for 
collisions of hydrogen atoms. Note that our predictions 
are not sensitive to a replacement of the atomic mass by 
the nuclear mass. 
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Formation of a tritium condensate in the high-field- 
seeking |0, 0) state would require some variety of non- 
magnetic trap, such as the dipole CO2 laser trap that 
has already produced a 87 Rb condensate M. Since the 
electric dipole polarizability of tritium is only 4.5 a.u., 
the 12 W C0 2 laser setup of § would only produce trap 
depths of the order of a few microkelvin. There appears 
to be no reason why much stronger CO2 lasers could not 
be utilized, however. Realistically, the formation of a 
tritium condensate will probably require laser intensities 
at least an order of magnitude more intense. The large 
magnitude of the zero-field scattering length would make 
evaporative cooling far more effective than is the case for 
spin-polarized hydrogen. Another possible way to cool 
the spinless substate of tritium would be to implement 
a recent proposal to cool an atomic gas through mag- 
netic field ramps across a Feshbach resonance 191. De- 
spite these technical difficulties to be overcome, recent 
improvements in trapping and cooling technology would 
appear to make the creation of a tritium condensate a 
viable possibility. 

To determine the bound state properties of (Tf )at dus- 
ters with N > 2, we first have to investigate the impor- 
tance of non-additive contributions to the many-body in- 
teraction potential. Although three-body contributions, 
i.e., the Axilrod- Teller term p0|| and three-body exchange 
terms are significant for the hydrogen trimer in its 
electronic ground state, they should be less important 
for the spin-polarized trimer since the classical atom- 
atom equilibrium distance of the quartet trimer is more 
than five times as large as that for the doublet ground 
state trimer. Our calculations show that inclusion of the 
damped Axilrod- Teller term |E| raises the ground state 
energy of the spin-polarized tritium trimer by roughly 
1.6%, and that of the larger clusters slightly more, e.g., 
by about 6 % for N = 40. To describe (Tf)jv clusters, 
we assume in the following a simple pairwise additive 
potential energy surface, V — J2i<j V( r ij)> where V(r) 
denotes the triplet 6 3 £+ two-body potential. Conceiv- 
ably, a more sophisticated many-body potential energy 
surface, which includes effects beyond the two-body po- 
tential, could modify our quantitative results somewhat, 
but we do not expect qualitative changes. 

For the Tf trimer, we use the adiabatic hyperspher- 
ical representation j22j] . Including only one adiabatic 
channel yields a single bound state with energy —1.60 x 
10~ 3 cm -1 . Coupling 25 adiabatic channels results in 
an energy of —3.19 x 10~ 3 cm -1 with an uncertainty of 
10~ 5 cm -1 and still no excited states. In the limit that an 
infinite number of channels are coupled, the bound state 
energy becomes exact, so the uncertainty is the result of 
including a finite number of channels. Upon inclusion 
of the Axilrod- Teller three-body term, the ground state 
energy is raised to —3.14 x 10~ 3 cm -1 — again with un- 
certainty in the last digit. There is thus a single L = 
bound state for quartet tritium (L is the total orbital an- 
gular momentum); no L > bound states are expected 



since none exist for the 4 He trimer |24| . 

Since the Tf dimer is unbound, the Tf trimer is a 
Borromean state (p5|. One may then ask: does the 
Tf trimer state have Efimov character ^(|? To inves- 
tigate this question, we apply a simple quantitative cri- 
terion (2^] , although others exist |28| . If the bound state 
disappears when one makes the potential more attrac- 
tive (here achieved by simply multiplying the two-body 
potential with an overall scaling factor greater than 1), 
then the state under investigation is an Efimov state; if 
the bound state does not disappear, it is not an Efimov 
state. Applying this criterion, our coupled-channel cal- 
culations indicate that the bound Tf trimer state is not 
an Efimov state. In short, the Tf trimer has exactly one 
L = bound state, a Borromean state that is highly 
diffuse spatially. 

To calculate the energetics and structural properties of 
(Tf )jv clusters with up to N — 40 atoms, we employ the 
diffusion quantum Monte Carlo (DMC) technique [p9| . 
This method solves the time-independent many-body 
Schrodinger equation essentially exactly, to within a sta- 
tistical error. Here, we employ the DMC method with 
importance sampling p9| , using a descendant weighting 
scheme p(| for the extrapolation of structural properties. 
Our guiding wave functions | p9[ , which enter the DMC 
calculation, have the analytical form given in Eq. (5) of 
Ref. and recover between 84 and 96% of the DMC 
energy when used in a variational quantum Monte Carlo 
calculation. The DMC technique as implemented here 
determines only the ground state of the system; accessing 
excited state properties is in general not possible without 
algorithmic modifications. 
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FIG. 3. Ground state energy per particle Eq/N for (Tf)jv 
clusters (pluses) and 4 Hejv clusters (diamonds) as a func- 
tion of N . Dotted lines are shown to guide the eye. In- 
set: Pair distribution P(r) for (Tf)s cluster (solid line) and 
4 He5 cluster (dotted line), calculated by the DMC technique 
with importance sampling. P{r) is normalized such that 
f™P(r)r 2 dr = 1. 
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For comparison, we find a DMC binding energy of 
—2.9(5) x 10 _3 cm _1 for the Tf trimer, in agreement with 
our hyperspherical calculation. The number in brackets 
denotes the statistical uncertainty. Figure || compares the 
ground state energy per particle Eq/N of (Tf)jv clusters 
(pluses) with those of 4 HeAr clusters (using the LM2M2 
potential jlj], diamonds). This shows that (TT)at clus- 
ters are even more weakly bound than 4 HeAr clusters with 
the same number of atoms. 

The comparison between (Tf )n and 4 HeAr clusters can 
be extended by considering their structural properties. 
For example, we find that the (Tf)s cluster has an av- 
erage interparticle distance of (r^-) of 22.2a.u.; the 4 Hes 
cluster, on the other hand, is (r^ ) = 13.6a.u. Even if one 
takes into account that the classical equilibrium distance 
of the tritium triplet potential is about 2.2a.u. larger 
than that of the He dimer potential, the difference be- 
tween the expectation values of the interparticle distance 
for these N = 5 clusters indicates that the spin-polarized 
tritium system is even more diffuse than the 4 Hes clus- 
ter. To illustrate this aspect further, the inset of Fig. || 
compares the pair distribution of the (Tf) 5 cluster (solid 
line) with that of the 4 He5 cluster (dotted line). Clearly, 
the pair distribution of the (T|)s cluster is much broader 
than that of the 4 He5 cluster. We find similar behavior 
for clusters with more particles. 

The tritium trimer is the smallest spin-polarized clus- 
ter (the dimer is unbound). As discussed above, the 
two-body potential for spin-polarized hydrogen is almost 
identical to that of tritium. The hydrogen atom, how- 
ever, is about a factor of three lighter. Consequently, 
there is no bound state for the H| trimer. An interest- 
ing question to ask is the following: how many atoms 
are needed to form a bound system of spin-polarized 
hydrogen atoms? Initial exploratory studies show that 
more than 100 atoms are needed to form a bound spin- 
polarized hydrogen cluster. The smallest spin-polarized 
hydrogen cluster could then be thought of as a "super- 
Borromean" cluster for which all smaller subsystems are 
unbound. A detailed study will be published elsewhere. 

In summary, this Letter proposes a number of intrigu- 
ing possibilities for the physics of tritium systems. We 
point out the possibility for forming a tritium condensate 
with controllable parameters via a Feshbach resonance. 
Down the road, one can imagine trapping an atomic hy- 
drogen gas together with an atomic tritium gas, or pos- 
sibly including deuterium to study fermion systems. We 
further found that the spin-polarized trimer possesses a 
Borromean, or halo, state. In addition, we mapped out 
the properties of larger "exotic" spin-polarized tritium 
clusters. Studies of spin-polarized clusters are interest- 
ing by themselves []32f| - as they enter new many-body 
physics regimes. For instance, unexpected physics may 
emerge from manipulating the two-body interaction via 
the Feshbach resonance in a cluster. In short, tritium of- 
fers a wealth of interesting physics by virtue of its weak 
attraction. 
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